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We study the quasiperiodic Harper's model in order to 
give further support for a possible universality of the criti- 
cal spectral statistics. At the mobility edge we numerically 
obtain a scale-invariant distribution of the bands S, which 
is closely described by a semi-Poisson P{S) — 4S'exp(— 25*) 
curve. The exp(— 25") tail appears when the mobility edge 
is approached from the metal while P{S) is asymptotically 
log-normal for the insulator. The multifractal critical density 
of states also leads to a sub-Poisson linear number variance 
T.2{E) oc 0.041_E. 
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During the past few years the statistical description 
of the energy levels in quantum systems has emerged 
as a new field of study, bringing together the areas of 
mesoscopic physics and quantum chaos [Q . In disordered 
metals, the states are extended with correlated energy 
levels characterized by level-repulsion found in appro- 
priate random matrix ensembles. This universal limit 
known as Wigner statistics also applies to quantum sys- 
tems with chaotic classical dynamics [EJ. In disordered 
insulators, the states are localized and the energy spec- 
tra obey Poisson statistics, which is the generic Hmit for 
intcgrable systems. At the metal-insulator transition, an- 
other scale invariant limit applies for the spectral fluctu- 
ations H] , associated with the multifractal nature of the 
critical eigenstates, which are neither extended nor local- 
ized H. However, the analysis is complicated by bound- 
ary effects g,g] which persist in the thermodynamic limit. 
But after taking specific boundary conditions, or aver- 
aging over a combination of them, the critical statistics 
resembles that obtained in certain weakly chaotic quan- 
tum systems, such as pseudointegrable rational triangle 
billiards M. Similar statistics also appears in certain 
many body spectra (e.g. two particles S in a disordered 
chain with on-site interactions), where P{S) is Wigner- 
like (linear) for small S and Poisson-like (exponential) 
at large S, overall described by the semi-Poisson curve 
P{S) = 4S'exp(— 25"). This result is accompanied by a 
sub-Poisson linear number variance 'S2{E) ~ x^ with 
X < 1. Semi-Poisson P{S) and sub-Poisson T^2{E) seem 
to be a common feature of different critical systems, rais- 
ing the question of a possible universality. The semi- 
Poisson P{S) can be exactly derived from a short range 
plasma model for the joint probability distribution of the 



energy levels, once the logarithmic pairwise interaction 
is truncated to nearest neighbours M. Screened inter- 
actions were also found to be consistent with the non 
universal long range part of the spectral fluctuations in 
disordered metals [p yic|] a nd characterize solvable models 
proposed by Gaudin |11| and Yukawa [Oj. Our aim is 
to show that the signatures of "critical chaos" appear at 
the mobility edge of the Id Harper model, supporting the 
idea of possible universality. 

The quasiperiodic Harper's equation |13|| has a mobil- 
ity edge in one dimension with no Wigner statistics for 
the metal, due to the quasi-ballistic nature of the ex- 
tended states. The Poisson law is also never reached for 
the insulator, due to the correlated localized states p4[ . 
At the critical point the spectrum is multifractal [[l5| with 
gaps G of all sizes distributed according to the inverse 
power-law P{G) ^ G'~(^+^' and D w 0.5, the spectral 
fractal dimension [ [l6|Jl7| . The distribution of gaps was 
interpreted as nearest-level spacing distribution charac- 
terized by strong level-clustering P{G ^ 0) ^ oo, which 
is very different from the level-repulsion ^(5* — > 0) ^ 
of the critical P{S) |||,|^. This apparently contradicts 
the idea of universality of the critical fluctuations and 
led us to revisit the issue. In order to obtain a possi- 
bly universal behaviour, we should (i) unfold the fractal 
spectrum to get rid of the average variation of the den- 
sity of states p{E) by keeping the rescaled local mean 
level spacing l/p{E) fixed, (ii) identify the crucial role 
of boundary conditions ^j^, and (iii) connect the crit- 
ical level statistics to the observed multifractality [[l8[ . 
The difficulty of (i) lies on the wildly fluctuating mul- 
tifractal spectrum and was circumvented in |lj] where 
p{E) was obtained by mapping the system onto a ring 
threaded by a flux and averaging over all possible fluxes 
(or equivalently by choosing different boundary condi- 
tions in the chain ends). This could suggest to average 
over boundary conditions in order to achieve universal- 
ity. Moreover, it leads us to argue that the distribution of 
bandwidths S is more appropriate than the distribution 
of gapwidths G for probing level statistics since we show 
that the normalized P{S) at the mobility edge is close to 
the scale-invariant semi-Poisson curve observed in other 
critical systems 0]. On the basis of these findings we 
suggest a certain universality of the critical fluctuations 
and reveal their connection to multifractality. 

We analyse the Harper's equation 



Elpn = XcOs{2'Kan + v)lpn + V'n+l + '^n-1, 



(1) 



where i/)„ indicates the particle wavefunction of energy 
E on site n, with i^ = ||l||. Eq. (1) for A = 2 also 
describes a particle in a two-dimensional lattice in a uni- 
form magnetic field with a flux quanta per elementary 
plaquette. When a is an irrational number the period 
of the potential is incommensurate with the lattice pe- 
riod. We can consider generic irrationals which cannot 
be approximated "too well" by rationals taking a as the 
limit of successive rationals M/N, so that the poten- 
tial becomes commensurate with the lattice with period 
N. Then the problem reduces to the diagonalization of 
N X N matrices with a spectrum which consist of N 
bands separated by A^ — 1 gaps. This also defines a scal- 
ing procedure where the incommensurate limit N —^ oo 
becomes equivalent with the thermodynamic limit. The 
states of Eq. (1) are extended when A < 2, with non- 
zero Lebesgue spectral measure 2|2 — A| independent of 
N. For A > 2 the states are localized and the spec- 
tral measure falls off exponentially fast. The most inter- 
esting case is the critical point A = 2 where the spec- 
trum and the wavefunctions are known to have hierar- 
chical multifractal structure [15 19-pl|]. Thouless showed 
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that the total band measure (suin of bandwidths) is pro- 
portional to 4. 6649744644. ..A^^^, which defines a uni- 
versal number p2]. The critical spectrum is multifrac- 
tal which implies that the various bandwidths S fall off 
with different power laws N~^/" , where a are singularity 
strengths having density f{a). At the edges of the spec- 
trum the bands scale as iV~^/"™'", the central band as 
7V-i/"™o. and the mean band is 4.6649744644.. .7V^i/"o, 
where a„im ~ 0.421, OLmax ~ 0.547 and ao = 0.5 are the 
most prominent spectral dimensions [ p5|p9| ]. 

We compute the bandwidths Si for long-A^ chains by 
finding the eigenvalues of the corresponding matrices 
for periodic and anti-periodic boundary conditions us- 
ing the symmetry of the potential A cos(27r crn) in order 
to diagonalize the symmetric and antisymmetric tridi- 
agonal matrices suggested in [g^. We find the non- 
overlapping bandwidths Si, i = 1,2, ..,7V and the gaps 
Gi, i = 1,2,..,A^— 1 for the inverse golden mean a — 
(v(5) — l)/2, which is the limit of the ratios of succes- 
sive Fibonacci numbers N . At the mobility edge A = 2 
the obtained normalized bandwidth distribution P{S) is 
displayed in Fig. 1 for various sizes N . It shows lin- 
ear level-repulsion and exponential tails, overall being 
described by the semi-Poisson P{S) = 4S'exp(— 5) curve, 
in sharp contrast with the power-law P{G) ^ G~^/^ ob- 
tained from the distribution of gaps |lq,0]. This sug- 
gests to examine whether the obtained P{S), which gives 
the distribution of the energy shifts when the boundary 
conditions are twisted for a commensurate potential of 
periodicity N, can be also understood as a meaningful 
spacing distribution for the incommensurate limit. 
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FIG. 1. (a) The normalized nearest level spacing distribution 
function at the mobility edge A = 2 for a characteristic size A'' = 
121393 (other sizes show very similar pictures). The solid curve 
is the semi-Poisson P{S) = 45exp(— 25). (b) P{S) in a semi-log 
plot for various sizes A^ with the semi-Poisson In P{S) = ln(45) — 
25" result. The computed total band multiplied by the given TV 
approaches the universal number of [fed] , within 4 to 8 decimals. 

The level fluctuations should be studied by consider- 
ing the distribution of (-E^+i — Ei)p{Ei), where p{E) is 
some averaged density of states at the energy E. In a 
disordered system this is usually done by averaging over 
the ensemble at E and in chaotic systems one can use a 
semiclasical formula for the mean density B. However, 
to unfold a multifractal spectrum seems a hopeless task 
since the level spacing fluctuates violently making diffi- 
cult to define a mean density. In other words one should 
consider the statistics of S/{G) around E, although the 
mean (G) cannot be defined. The solution was suggested 
in rcf. |1J] which uses the unfolded levels from a mean 
p{E) (or (G)) obtained by repeating periodically the iV- 
site chain or changing the flux 0, which describes the 
Bloch boundary condition 4'n+N — exp(z0)^„. A sim- 
pler way to do this is by noticing that the number of 
levels filling every band (which can be due to repeating 
periodically the long chain of N sites to infinity [Q), is 



fixed {^ 1/-^) since to a given (j) correspond N levels, 
one in every band. Moreover, each band has a smooth 
density of states which can be made constant indepen- 
dently of its form, by scaling the bandwidths Si with the 
fixed N. Thus, every bandwidth contributes a delta func- 
tion to P{S) centered at Si, similarly to a delta function 
P{S) at the mean spacing obtained for a smooth peri- 
odic band and the level-spacing statistics can be associ- 
ated with the bandwidth distribution, essentially without 
unfolding. Viewed this way P(5') concerns the distribu- 
tion of the mean spacings, each one of them arising from 
the levels filling every band. The shown scale-invariant 
P{S) described by the semi-Poisson curve (Fig. 1) can 
also justify the results at the metal-insulator transition 
in three dimensional disordered systems, where a similar 
P{S) was obtained by averaging over boundary condi- 
tions Ig,^. 

A remarkable connection was recently made between 
level-statistics and the multifractality of critical wave- 
functions ||l^. At the mobility edge the long range 
spectral fluctuations are described by a linear number 
variance T,2{E) = xo + X^^ with level compressibility 
X = {1/2){1 — D2 /d) related to the correlation dimension 
D2 of the critical wavefunctions and the space dimension 
d. The quasiperiodic model has no unique wavefunction 
exponent D2 to describe the compressibility, since Dj 
varies with the critical energies. In the Harper's model 
T,2{E) is directly influenced by the spectral fluctuations 
summarized in the presence of the multifractal dimen- 
sions Dq,q G (—00,-1-00) and their Legendre transform 
a— f{a) ||l^,|l^. They are locally defined around the 
spectral energy Ein, via 



M{E,n+E)-U{E,n)^E° 



(2) 



where Af{E) is the number of levels from the lowest en- 
ergy up to E. 



In the formalism of multifractals a spectral scaling 
function r(q) is usually evaluated by box counting in en- 
ergy, where an exact partioning of boxes is suggested by 
the bandwidths Si since no levels fall in the gaps. The 
scaling dimensions are obtained by scaling the moments 
of the computed distribution of Si [O via 



N 



E^^ 



Mq) 



N". 



(3) 



For up to A^ = 313897 it gives the spectral exponents 
Dg = T{q)/{q - 1) with D-i = 0.5, L>o = £> « 0.498, 
Di sa 0.496, D2 ~ 0.493, etc. Their Legendre transform 
gives the singularities a = — ^^ and their density f{a) — 
qa — T{q) with the extreme dimensions -D+oo = ctrnim 
£'_oo = <^max- The position of the maximum lies at 
ao = 0.5 with f{ao) — D (see inset of Fig. 2). 
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FIG. 3. The nearest level spacing distribution function in the 
metallic phase A < 2 for N = 121393. The continuous line is the 
semi-Poisson limit approached when A = 2. 

The level number variance 
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S2(£;) = / dEU^{E,n +E)- M{E,n)f 



dE,^{U{E„,+E)-U{E„,))] , (4) 



is defined by taking energy boxes of size E over all pos- 
sible Ein in the spectrum. In the box multifractal for- 
malism the knowledge of a — f{a) makes Eq. (4) strictly 
equivalent to 



S2(S) = 



daf{a)E^°' - 



daf{a)E" 



(5) 



FIG. 2. The number variance at the mobility edge A = 2. Inset: 
the a-f{a) spectral dimensions. 



which is approximately linear E2 ~ 0.012 -I- 0.041_E as 
shown in Fig. 2, leading to level-compressibility x ~ 



0.041. The obtained x should be contrasted with the 
higher values x ~ 0.10,0.125 for d = 2 and x ~ 0-32 
for d = 3 critical disordered systems [Q. In the one- 
dimensional mobility edge x is smaller, inversely propor- 
tional to d as expected from the formula for x via the 
wavefunction exponent D^ Pq]. The computed value of 
X applied to this formula gives ZJj ~ 0.918, which can 
be compared with Dj ~ 0.82 obtained for states near 
the band center [£lj . Since all the levels in the band cor- 
respond to critical states with variable Z?2 we have also 
computed directly the average exponent for all critical 
states which gave the much lower value Z?2 ~ 0.537. 

The level statistics of the Harper's equation for A < 2 
or A > 2 can be simply understood from a duality ar- 
gument [Q. In the metallic regime we follow P{S) 
for increasing A in Fig. 3. The Poisson tail exp(— 25*) 
is seen to emerge when approaching the critical point 
A = 2. The spectral dimensions in the metallic limit 
are simply Dq = 1, for q < 2, and Dq = 1/(2 — 2/q), 
for q > 2. For example, for A = we trivially obtain 

N, which leads to 
3, with only one 



{2n/N)\sin{'KJM/N)\,j = 1,2,, 
shown in Fig. 



singularity a — 0.5, f{a ~ 0.5) = 0. More and more 
singularities develop as we approach the critical point 
A = 2 where a lies in {amin,cimax) with finite density 
/(a). The obtained critical semi-Poisson P{S) gives an- 
other hint towards the fascinating problem of the analyt- 
ical computation of the spectral multifractal dimensions, 
possibly from finite- iV corrections to strings in the recent 
formulation via Bethe-anzatz equations |23| . In the local- 
ized regime A > 2 the bands become exponentially small 
and P{S) asymptotically log-normal since the localized 
levels are spatially correlated. In the extended phase al- 
most all a = 1 and in the localized phase a = 0, which 
lead to trivial T,2{E) = in the two limits. 

The main result of our study, which is able to treat nu- 
merically very long chains, is the appearance of a semi- 
Poisson P{S) = 4S'exp(— 25") distribution at the mobil- 
ity edge in one dimension. This associates the mobil- 
ity edge with level-repulsion for small spacings, rather 
than level-clustering, in agreement with 3d critical disor- 
dered systems. At the mobility edge both the multifrac- 
tal spectrum and the nearest-level statistics are seen to 
arise from the semi-Poisson distributed bandwidths Si, 
whose moment scaling describes a fractal spectrum with 
infinite many singularities a — f{a). A second result di- 
rectly obtained from the multifractal spectrum is a sub- 
Poisson linear number variance Y,2{E) ~ xE^ in addition 
to the approximate validity of the proposed formula for 
X in terms of the multifractal wavefunctions. In con- 
clusion, the obtained universal critical statistics for the 
Harper's model verifies numerically the "critical chaos" 
scenario summarized in the semi-Poisson P{S) and the 
sub-Poisson linear number variance T,2{E). "Critical 
chaos" is seen to be intimately connected to spectral and 



wavefunction multifractality. Our study can also shed 
light on the average over boundary conditions in the crit- 
ical statistics of disordered systems. 
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